On the compressibility of deformable spin chains in a vicinity of quantum critical 

points 



Oleg Derzhko, 1,2,3 Jozef Strecka, 4 and Lucia Galisova 5 

'institute for Condensed Matter Physics, National Academy of 
Sciences of Ukraine, 1 Svientsitskii Street, L'viv-11, 79011, Ukraine 
2 Department for Theoretical Physics, Ivan Franko National University of L'viv, 12 Drahomanov Street, L'viv-5, 79005, Ukraine 
s Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 1-34151 Trieste, Italy 
4 Department of Theoretical Physics and Astrophysics, Faculty of Science, 
P. J. Safdrik University, Park Angelinum 9, 040 01 Kosice, Slovak Republic 
5 Department of Applied Mathematics and Informatics, Faculty of Mechanical Engineering, 
Technical University, Letnd 9, 042 00 Kosice, Slovak Republic 
(Dated: October 16, 2012) 

We calculate the ground-state compressibility of a deformable spin- 1/2 Heisenberg-Ising chain 
with Dzyaloshinskii-Moriya interaction to discuss how a quantum critical point inherent in this 
spin system may manifest itself in the elastic properties of the underlying lattice. We compare 
these results with the corresponding ones for the spin- 1/2 Ising chain in a longitudinal or transverse 
field and the spin- 1/2 XX chain in a transverse field. The inverse compressibility of the spin- 1/2 
XX chain in a transverse field exhibits a hysteresis in a vicinity of quantum critical point that is 
accompanied with the finite jump of compressibility. Contrary to this, the inverse compressibility 
diminishes continuously close to a quantum critical point of the spin-1/2 Ising chain in a transverse 
field and the spin-1/2 Heisenberg-Ising bond alternating chain. 
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I. INTRODUCTORY REMARKS 

Quantum spin chains provide an excellent play-ground 
for theoretical studies of collective quantum phenom- 
ena. They may exhibit very rich low-temperature be- 
havior showing numerous quantum phases (sometimes 
exotic ones which are characterized by non-conventional 
order parameters) and quantum critical points^. More- 
over, including spin-lattice interactions may lead to new 
phenomena like, e.g., spin-Peierls instabilities^. Impor- 
tantly, quantum spin chains in many cases are amenable 
to rigorous analysis that permits to obtain results free of 
using any uncontrolled approximations^—. At the same 
time, real solid-state representatives of exactly solved 
quantum spin-chain models have become available due 
to a recent progress in material sciences and comparison 
of theoretical predictions with experimental observations 
has thus become possible^. 

In the present paper we wish to discuss how quantum 
critical point, inherent in a system of interacting quan- 
tum spins, may manifest itself in the properties of the 
underlying lattice. For this purpose we consider sim- 
ple but nontrivial exactly solvable spin models focusing 
in particular on a spin-1/2 Heisenberg-Ising chain with 
Dzyaloshinskii-Moriya interaction. We also consider a 
harmonic one-dimensional lattice and use an adiabatic 
treatment assuming that the lattice may rearrange itself 
to suit a "request" of the interacting quantum spins. The 
results to be obtained will be confronted with the cor- 
responding ones for three exactly soluble quantum spin 
chains, more specifically, the spin-1/2 Ising chain in a lon- 
gitudinal or transverse field and the spin-1/2 XX chain 



in a transverse field. In our study we calculate the lat- 
tice compressibility to show how its behavior indicates 
peculiarities in those quantum spin systems. Our goal is 
not to describe any specific material, but to analyze ac- 
curately model systems amenable to exact calculations. 
We believe that our analysis of correlation between lat- 
tice properties and peculiar behavior of spin system will 
have some qualitative merit in general. 

The paper is organized as follows. First, in Sec. |TT] 
we illustrate a calculation scheme considering some well- 
known examples of exactly solved quantum spin chains 
including the spin-1/2 Ising chain in a longitudinal field, 
as well as, the spin-1/2 XX chain in a transverse field. 
In Sec. IIIII we examine the elastic properties of the 
deformable spin-1/2 Ising chain in a transverse field. 
Our main emphasis is laid on the deformable spin-1/2 
Heisenberg-Ising chain with Dzyaloshinskii-Moriya inter- 
action considered in Sec. IIVI In this section we exam- 
ine the compressibility of the spin lattice by discussing 
its relation to a quantum critical point inherent in the 
spin system. Finally, our findings will be summarized in 
Sec. ED 



II. COMPRESSIBILITY OF SOME QUANTUM 
SPIN-CHAIN MODELS 

In this section, we discuss a calculation of the com- 
pressibility of a lattice with interacting spins. First, con- 
sider an "empty" (i.e., without spins) one-dimensional 
lattice of N sites which is perfectly periodic and has 
minimum energy for some lattice constant do- Apply- 
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ing pressure along the lattice we may either elongate or 
shorten it. We assume that the extension or reduction 
of the lattice length is uniform, i.e., the change of dis- 
tance between any two neighboring sites is in equilib- 
rium a — ao oc S, and the total elastic energy penalty is 
NaS 2 /2, where a is the (bare) elastic constant of the lat- 
tice. We consider the variational Gibbs free energy per 
site g(T,p; 5), or more precisely, the variational enthalpy 
per site g(T = 0, p; S) since we restrict ourselves in what 
follows to the zero-temperature case 

g(T = 0,p;5) = ^aS 2 +pS, (2.1) 

where p is the (dimensionless) pressure. The equilibrium 
value of S denoted further as S(p) is obtained by mini- 
mizing the variational enthalpy g(T = 0,p;<5) (|2 . 1 [) with 
respect to S (see Refs. llollTTh . i.e., by solving the equation 
= dg(T = 0,p;5)/dS = a8 +p, with the result 

*(P) = ~. (2-2) 
a 

Note that S(p) > (S(p) < 0) [i.e., stretching (shrinking) 
of the lattice occurs] if p < (p > 0). Substituting 
5 = S(p) (|2~2|) into Eq. (|2J|) we obtain the enthalpy 
per site g(T = 0,p) = g(T = 0,p;S(p)) = -p 2 /(2a). 
Moreover, using Eq. ()2.2|) we can calculate the observed 
lattice compressibility (or the inverse elastic constant) 

i - - (m = i. ( 2 .3) 

x \ op J a 

Obviously, we have found that the observed elastic con- 
stant of the lattice x coincides with the bare elastic con- 
stant a as it should for the case at hand. 

Now we populate this lattice with N interacting quan- 
tum spins s = 1/2 which are represented by Pauli matri- 
ces. We assume only the nearest-neighbor exchange in- 
teractions between the spins which obviously depend on 
the intersite distance. Thus, Eq. (|2.1[) will now contain 
on the right-hand side a <5-dependent term eo(d) corre- 
sponding to the ground-state energy of the spin system 
g(T = 0,p; S) = aS 2 /2 + pd + eo(S), and we may argue 
about the effect of the spin system on the observed lattice 
compressibility 1/x defined in Eq. (|2.3[) . It is worthwhile 
to consider briefly several well-known spin models. 

Ising chain in a longitudinal fields. We begin with 
the spin Hamiltonian 

H = J2 ( J «+i ~ K) > h>0 (2.4) 

n 

and assume J = Jq(1 — kS) with k = 1 for simplicity. 
The adopted linear dependence of J on S can be justified 
only if S < 1. 

For the ferromagnetic exchange interaction Jo = 
— | Jo | < the ground-state energy per site is given by 
eo(S) = — | Jo | (1 — <5)/4 — h/2 that results in a new equi- 
librium value of 5, 8(p) = — (p + | Jo|/4)/a, but leaves 



the inverse compressibility unchanged, x = a. {Note 
that for p = we have 5(0) = — |Jo|/(4a), i.e., the 
lattice shrinks and the ferromagnetic exchange constant 
J = Jo[l + | Jo|/(4a)] becomes larger than J .} 

For the antiferromagnetic exchange interaction Jq = 
| Jo | > 0, the ground-state energy per site is given cither 
by e 0w (<5) = -I Jo|(l - *)/4 if h< \ J |(l - S) (weak-field 
regime) orbye 0s ((5) = \ J \(l-S)/i-h/2i£h> \Jo\{lS) 
(strong- field regime). Let us introduce the critical fields 
h c /\Jo\ = 1 +p/a, h cl /\J \ = h c /\J \ - \ J \/(4a), and 
h C 2/\Jo\ = h c /\ Jo | + 1 Jo|/(4a), the critical pressures p c = 
(h/|Jo|-l)a, Pd = p c -\ Jo|/4, andp c2 = p c + \ Jo|/4, and 
the two equilibrium values of S, 5\(p) = — (p + \ Jo|/4)/a 
and 5 2 {p) = -(p - | Jo |/4) /a. 

Considering only the global minimum of the variational 
enthalpy g(T = 0,p,h;S) one arrives at the following 
conclusions, i) At fixed values of p, the equilibrium value 
of 8 is either <5i(p) for h < h c or S^ip) for h > h c and 
there is a jump in between these two equilibrium values 
at the critical field h = h c . ii) Under the constant field 
h, the equilibrium value of S is ^(p) for p < p c , Si(p) 
for p > p c and the abrupt change in S (p) appears at the 
critical pressure p = p c . 

However, it should be also stressed that the depen- 
dence of variational enthalpy g(T = 0,p,h;S) on 8 may 
exhibit an additional local minimum besides the global 
one when driving the antiferromagnetic Ising chain suf- 
ficiently close to the critical field h c or the critical pres- 
sure p c . This implies a possibility of observing hysteresis 
phenomena originating from the first-order phase tran- 
sitions driven either by varying of field h or pressure p. 
Considering the metastable states, which correspond to 
the local minima of the variational enthalpy, one arrives 
after a simple analysis at the following conclusions, i) 
For fixed p, the equilibrium value of 5 is 6\ (p) when h 
increases from to h c2 and then it jumps to S 2 (p) when 
h further increases (increasing-field regime). However, 
the equilibrium value of 8 is S 2 (p) when h decreases from 
oo to h c \ and then it jumps to 5\(p) upon further de- 
crease of h (decreasing-ficld regime). Furthermore, the 
relations d5\(p)/dp = dS 2 (p)/dp = — 1/a give x = a. ii) 
For fixed h, the equilibrium value of S is 5 2 (p) when p 
varies from — oo to p 2c and then it jumps to 5i(p) as p 
exceeds p 2c (increasing-pressure regime). On the other 
hand, the equilibrium value of 6 is 8\(p) when p varies 
from oo to p\ c and then it jumps to 6 2 (p) whenever p 
becomes less than p\ c (decreasing-pressure regime) . Fur- 
thermore, x — a unless p ^ p c2 but x(p c2 ) = for the 
increasing-pressure regime, whereas x = a unless p =/= p c \ 
but x(p c \) = for the decreasing-pressure regime. 

XX chain in a transverse field. We pass to the one- 
dimensional spin-1/2 system with isotropic XY interac- 
tion in the presence of a transverse magnetic field. The 
Hamiltonian of the spin model reads 

71 

and J = Jo(l — S). The sign of J is not important for 
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the thermodynamic quantities and therefore, we may set 
J > without loss of generality. Although this model 
in t he p resent context was discussed in some detail in 
Ref. [ll], we report here a brief account of to some extent 
known results for self-consistent readability of our paper 
and for reader's convenience. The magnetic contribution 
eo{5) to the variational enthalpy can easily be found using 
the Jordan- Wigner fermionization approach^. Separat- 
ing the weak-field regime and the strong-field regime we 
have 



regime). The inverse compressibility for the increasing- 
field regime follows 



e (S) 
eow(<5) 
eo s (£) 



e 0w (8), h < J (l - 8), 
e 0s (5), h>J (l-5), 

Mi - s) 



sin kv, 



/; 
2 



(2.6) 



with cos/cf = h/[J (l — 8)]. Again, g(T = 0,p 7 h;S) 
versus 8 dependence can have two minima suggesting a 
possible existence of first-order phase transitions accom- 
panied by hysteresis phenomena. 

More specifically, consider the equation for S(p) in the 
weak-field regime 



aS + p + 



Jo 



h 2 



J 2 (l-<5) 2 



(2.7) 



Fixing p and varying h from to oo one finds that i) 
Eq. ()2.7[) has only one solution which corresponds to a 
minimum of g w (T = 0,p,h;8) = aS 2 /2 + pS + eo w (S) 
for any h < h c \ = (1 + p/a)Jo; ii) at h = h c \ a sec- 
ond solution (at —p/a) corresponding to a maximum of 
<7w(T = 0,p,h;8) emerges; Hi) Eq. (|2.7j) has two solu- 
tions for h c i < h < h C 2 that are gradually approaching 
each other as h increases; iv) at h = h C 2 two solutions 
coincide and Eq. (|2.7j) has no appropriate solutions for 
larger h. Let us denote the smaller root of Eq. (|2.7[) by 
5i(p). The equation for d(p) in the strong- field regime 
yields 82 ip) = —p/cx which corresponds to a minimum of 
g s (T = 0,p,h;6) = a5' 2 /2 + p5 + eo s (S). One more char- 
acteristic field, h c i < h c < h c2 , is determined from the 
condition g w (T = 0,p, h c ; Si(p)) = g s (T = 0,p, h c ; S 2 (p)). 

Now one may immediately conclude that the equi- 
librium value of 8 is Si(p) for h < h c and 82 ip) for 
h > h c when considering the thermodynamically sta- 
ble solution corresponding to the global minimum of the 
variational enthalpy only. On the other hand, the inclu- 
sion of metastable states into our analysis allows one to 
study an interesting hysteresis phenomena resulting from 
the magneto-elastic coupling. Under this circumstance, 
one actually finds that the equilibrium value of 8 is 8\ip) 
when h increases from to h C 2 and then it jumps to 
8zip) > &i{p) when h further increases (increasing- field 
regime). However, the equilibrium value of <5 is 62 ip) 
when h decreases from 00 to h c \ and then it jumps to 
Si(p) < $2(75) upon further decrease of h (decreasing- field 



x = a — 



tt (1 - <mp)) v^a-^b)) 2 -/* 2 
h 2 

~ a+ tt 2 (aS 1 (p)+p)(l-S 1 (p)) 3 



(2.8) 



until h < h C 2, but xr = a if h > h C 2- For the decreasing- 
field regime x = a until h > h c \ but n follows Eq. (|2.8[) 
if h < h cl . 

Similarly, fixing h and varying p from —00 to 00 one 
finds that i) Eq. (|2.7j) has no appropriate solutions for 
p < p c i] ii) Eq. (|2.7p has two solutions if p exceeds 
Pd and the smaller one, 5\{p) 1 corresponds to a min- 
imum of g w (T = 0,p,h;5); Hi) the larger solution (at 
—p/a) corresponding to a maximum of g w (T = 0,p, h; 5) 
disappears at p = p C 2] iv) Eq. (|2.7[) has only one so- 
lution <5i(p) for p > p C 2- One more characteristic pres- 
sure, pd < p c < p C 2, is determined from the condition 
g w (T = 0,p c , h; Si(p c )) = g s (T = 0,p c , h; S 2 (p c ))- 

Considering only the stable solution with the low- 
est variational enthalpy one finds that the equilibrium 
value of 5 is either 8 zip) if p < p c or 8\(p) if p > p c . 
It is also quite obvious from the aforementioned anal- 
ysis that the striking pressure-induced hysteresis phe- 
nomenon emerges when taking into consideration an ex- 
istence of the metastable states. Indeed, it can be readily 
proved that the equilibrium value of 8 is 82 ip) when p in- 
creases from to p C 2 and then it jumps to 8\ip) < 82 ip) 
when p further increases (increasing-pressure regime). 
On the other hand, the equilibrium value of 8 is given by 
8\ip) when p decreases from 00 to p c \ and then it jumps 
to <!>2(p) > 8\ip) upon further decrease of p (decrcasing- 
pressure regime). Furthermore, x — a for the increasing- 
pressure regime until p < p C 2 but k follows Eq. (|2.8[) if 
P > Pc2- For the decreasing- field regime the inverse com- 
pressibility follows Eq. (|2.8|) if p > p c \ but x = a if 

P < Pel- 

In Figs. [T] and [5] we demonstrate dependences of the 
inverse compressibility x on h and p along with corre- 
sponding dependences of the equilibrium value of 8 de- 
picted in the insets. It is noteworthy that the displayed 
figures may serve in evidence of the aforedescribed ficld- 
and pressure-induced hysteresis phenomena and besides, 
one also observes here an abrupt jump of x and Sip) at 
the first-order phase transitions. The interested reader is 
referred to the paper by Orignac and Citroii for a more 
complete discussion of the dcformable spin- 1/2 XX chain 
in a transverse field. 

A brief summary of this section is as follows. Although 
the equilibrium value of 8ip) may vary after placing in- 
teracting spins on a lattice, the observed elastic constant 
x does not necessarily feel the appearance of interact- 
ing spins as for instance in the spin-1/2 ferromagnetic 
Ising chain in a longitudinal field. Contrary to this, the 
spin-1/2 XX chain in a transverse field clearly shows a 
ground-state peculiarity of the quantum spin model in 
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FIG. 1: (Color online). Inverse compressibility x versus field 
h for the spin-1/2 XX chain in a transverse field (|2.5[) with 
Jo = 1, a. = 2 at p = 0. The line with up-pointing arrow (in 
red) corresponds to the increasing-field regime, the line with 
down-pointing arrow (in blue) corresponds to the decreasing- 
field regime, and the line without an arrow (in green) cor- 
responds to the global minimum of the variational enthalpy. 
Inset: Corresponding S(p = 0) versus h dependences. 
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FIG. 2: (Color online). Inverse compressibility x versus pres- 
sure p for the spin-1/2 XX chain in a transverse field (|2.5|) 
with Jo = 1, a = 2 at h = 1.1. The line with down-pointing 
arrow (in red) corresponds to the increasing-pressure regime, 
the line with up-pointing arrow (in blue) corresponds to the 
decreasing-pressure regime, and the line without an arrow (in 
green) corresponds to the global minimum of the variational 
enthalpy. Inset: Corresponding S(p) versus p dependences. 



the elastic properties of the underlying lattice: the in- 
verse compressibility x noticeably diminishes while the 
spin chain approaches a quantum critical point indicat- 
ing a change from the spin-liquid phase to the ferromag- 
netic phase. More precisely, the inverse compressibility x 
shows a jump to a finite value accompanied with hystere- 
sis. Formally, a specific behavior of the elastic properties 
arises from a specific dependence eo [S) for the spin model 
at hand at a quantum critical point. 



III. ISING CHAIN IN A TRANSVERSE FIELD 

In this section, we consider the one-dimensional spin- 
1/2 transverse Ising chain given by the Hamiltonian 



H 



E 



hs z n ) , h > 



(3.1) 



with J = Jo (1 — 8). Again we may assume Jo > without 
loss of generality. This model, although to our knowledge 
has not been considered so far, can be again studied rig- 
orously using the Jordan- Wigner fermionization^^. The 
magnetic ground-state energy per site reads 



e (S) = - 



J (l -S) + 2h 
2^ 
Jo(l - S) 



E(z), 



2h 



n 2 



J (l -S) + 2h 



(3.2) 



where E(z) = J^ 2 d(j>\/l — z 2 sin 2 <j> is the complete el- 
liptic integral of the second kind with the modulus z de- 



fined through Eq. (|3.2[) . Recal l 13 ' 14 that around z sa 1 we 
haveE(z) « 1 + (l/2)[ln(4/Vl - z 2 ) - 1/2](1 - z 2 ), that 
implies possible peculiarities when Jo(l — 5) — 2h = 0. It 
can be readily verified that the dependence of variational 
enthalpy g(T = 0,p,h;8) = a8 2 /2 + p8 + e (S) on the 
distortion parameter 5 has only one minimum at 5(p), 
whereas the value of 8{p) obeys the condition 



aS + p - 



de (S) 
88 



0, 



and the inverse compressibility then follows from 
a 2 e (<5) 



x = a + 



d8 2 



(3.3) 



(3.4) 



5=6(p) 



Here, eo(8) is given by Eq. (|3.2[) . 

In Figs. [3] and U] we demonstrate dependences x(h) at 
fixed p and x(p) at fixed h along with corresponding de- 
pendences of the equilibrium value of 8 as they follow 
from the numerical solution of Eqs. (|3.3p and (|3.4p with 
eo(S) (|3.2p . which have been solved by subsequent numer- 
ical integration and differentiation. Clearly, if the mini- 
mum of g(T = 0,p,h]8), which occurs at 8(p) given by 
Eq. (|3.3p . satisfies the relation Jq(1 — S(p)) —2h = 0, i.e., 
when 5(p) = 8 C with 8 C = 1 — 2h/ Jq, peculiarities in vari- 
ous quantities may be expected. For fixed p and varying h 
we find the critical field h c /J = l/2+\p+J /{2n)]/(2a) 
at which S(p) =S C = l-2h c /J = -\p+J /(2ir)]/a. For 
fixed h and varying p we find the critical pressure p c = 
-a(l-2/i/J )-Jo/(27r) at which 5(p c ) = 8 C = l-2h/J . 
It is noteworthy that the most obvious changes in the 
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FIG. 3: (Color online). Inverse compressibility yt versus field 
h for the spin-1/2 Ising chain in a transverse field (|3.ip with 
Jo = 1, a = 2 at p = —0.1 (green line), p = (red line), 
and p = 0.1 (blue line). Inset: Corresponding 8(p) versus h 
dependences. 




FIG. 4: (Color online). Inverse compressibility ^versus pres- 
sure p for the spin-1/2 Ising chain in a transverse field (|3.1[) 
with Jo = 1, a — 2 at h = 0.47 (green line), h — 0.5 (red line), 
and h = 0.53 (blue line). Inset: Corresponding 5(p) versus p 
dependences. 



monotonous 8{p) dependence on h (on p) at fixed p (at 
fixed h) appear close to the critical field h c (critical pres- 
sure p c ), see the insets of Figs. [3] and |3J 

Let us take a closer look on Eqs. (|3.2[) - (|3.4[) . For this 
purpose we rewrite Eq. (|3.3p in the form 



aS + p 

where K(z) = J^ 2 dtp/ \J\ — z 2 sin 2 <f> is the complete 



elliptic integral of the first kind which behaves as K(z) 
ln(4/\/l — z 2 ) around z 



l-I&M, and Eq. (031) in the 



form 



J (l-5) + 2ftdE(«) 



2?r(l - J) 2 
J (l -S)-2h dK(z) 



d: 



4tt(1 - 6) dz 



4tt(1 - 6) 



dz 



3.6) 



d~E(z)/dz = [E(z) - 
■E(z)/[z(l - z 2 )], and 



Here, z is given by Eq. (|3.2j) 
K(z)]/z, dK(z)/dz = -K(z)/z- 
dz/dS = 4:J h[J (l~6)-2h}/{z[J (l-S)+2h} 3 }. Clearly, 
Eqs. (|3.5[) and p.6|) immediately yield the dependences 
p(8) and x(5), respectively. 

Consider now, for instance, the case when h is fixed 
and p varies. Inverting the dependence p = p(5) which 
follows from Eq. (|3.5|) we may obtain S(p) (see the inset 
in Fig. |4}, whereas combining the dependences p = p(5) 
and yc = x(S) which follow from Eqs. (|3.5j) and (|3.6[) we 
may obtain x{p) (see the main panel in Fig- HJ) - 

If 6 = S c + e is close to 5 C (i.e., |e| — > 0), \J\ — z 2 sa 
\e\J /(4:h) and Eq. (|3.5|l becomes 



.In 



16h 



- 1 



8tt/i 



T2 

8tt/i 



In lei 



-C-^-eln|e| 



InC 



8nah 
•'u 



hi 



I6/1 



1. 



(3.7) 



Although |e| in Eq. f|3 . T[) is small, it is not less than 
C (for ft/ J = 0.45... 0.55, a/J = 2 we have C w 
3.98 • 10- 10 ...3.19 • 10~ 12 ), i.e., |e| > 1, in order not 
to violate the relation sgn(p — p c ) = — sgn(e). Differ- 
entiating Eq. (|3.7p with respect to 5 we obtain k = 
[Jg/(87r/i)](ln |e| + 1) and yt falls down to J%/{8nh) at 
|e| = 1. However, these values of k cannot be reached 
numerically, cf. Fig. 01 

A brief summary of this section is as follows. Both 
the equilibrium value of S(p) and the observed elastic 
constant yc clearly show a ground-state peculiarity of the 
quantum spin model put on the lattice. In particular, the 
inverse compressibility yt is noticeably diminished while 
the spin chain approaches a quantum critical point in- 
dicating a change from the Ising phase to the paramag- 
netic phase. In contrast to the transverse XX chain, the 
inverse compressibility yc for the transverse Ising chain 
decreases continuously but not discontinuously. 

Now we apply the elaborated scheme to the spin-1/2 
Heisenberg-Ising bond alternating chain refined by an ad- 
ditional Dzyaloshinskii-Moriya interaction. The model 
exhibits a quantum phase transition point driven by a 
relation between the interaction constants, which sepa- 
rates the disordered phase (weak Ising interaction) and 
the long-range ordered phase (strong Ising interaction), 
see Rcf.|5|. Our task is to follow how this quantum critical 
point manifests itself in the elastic properties. 
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IV. HEISENBERG-ISING CHAIN WITH 
DZYALOSHINSKII-MORIYA INTERACTION 

In this section, we focus on the ground-state proper- 
ties of a regularly alternating spin- 1/2 antiferromagnetic 
Heisenbcrg-Ising chain with a Dzyaloshinskii-Moriya in- 
teraction. The chain consists of N — 27V sites with spins 
s = 1/2. Moreover, the XX Z Heisenberg interaction 
and the z-component of the Dzyaloshinskii-Moriya inter- 
action regularly interchange with the Ising interaction in 
this chain. The Hamiltonian of the spin model reads 



H 



E 

m—l 



H 2 



with 



H 



2m— 1,2m 



D 



2m -I s 2m 



'2ro— l°2m 



m,2m+l J 



S L-1 S L + ^ S 



'2m-r2raJ 



(4.1) 



2m— l°2m) 



and 



H 



2m,2m+l — 3\ s 2m s 2m+\- 



(4.2) 



(4.3) 



Here the parameters Jh > and D denote the anti- 
ferromagnetic Heisenberg and the Dzyaloshinskii-Moriya 
interactions between 2m— 1 and 2m spins, A controls the 
anisotropy of the Heisenberg interaction, and the param- 
eter Ji > denotes the antiferromagnetic Ising interac- 
tion between 2m and 2m + 1 spins. We impose periodic 
boundary conditions for convenience. 

Although a merit of the introduced quantum spin chain 
is possibility to perform accurate calculations, it might be 
worthy to notice here that mixed Heisenberg-Ising inter- 
actions do occur in some magnetic compounds containing 
lanthanidc ionsi^. Moreover, the magnetic subsystem of 
such compounds may have one-dimensional geometry. 

The Hamiltonian (|4.1j) . (|4.2p . (|4.3j) corresponds to 
the exactly solvable Heisenberg-Ising model proposed in 
Ref. i (see also Refs. EMI). The calculation of the 
ground-state energy of the model in brief looks as fol- 
lows. First we eliminate the Dzyaloshinskii-Moriya term 
from the Hamiltonian after performing an appropriate 
spin-coordinate transformation^—, 



s 2m ~> s 2m COS (f + s| m sin tf, 
S 2m -> ~ S 2m sin <P + 4m C0S <P 



(4.4) 



with tan if = D/ Jh for m = 1, . . . ,7V, which results 

tl2m-\,2m = J XY \ s 2m-l s 2m + s 2m-\ S 2m) 



J H As 



2m-l 2m! 



JxY — 



D 2 . 



(4.5) 



From Eq. (|4.5[) we see that the introduced 
Dzyaloshinskii-Moriya interaction effectively increases 
the XY component of the Heisenberg interaction with 
respect to the Ising component of this interaction. 



Then we follow the arguments of Refs. mm The 
ground state of the antiferromagnetic model (|4.1[) , (|4.5p . 
(|4.3[) must lie in a particular subspace, where all Heisen- 
berg bonds are in one of two states: either (| J,)2m-i| t 

>2m + | t)2m-l| |>2m)/V2or (| |) 2 m-l| t)2m~| t>2m-l| I 

)2m)/v / 2- Introducing the raising and lowering operators 
for the mth pair, a m and a m , m = 1, . . . ,7V, which satisfy 
the Fermi commutation relations if they are attached to 
the same site, but obviously commute if they are attached 
to different sites, the Hamiltonian (|4.1j) , (|4.5|) . (|4.3[) be- 



comes 



^7 ^7 A ^ 

m—l 



"J Wm + a m) (4i+l + 



(4.6) 



After applying the Jordan- Wigner transformation^, 

ci = oi, c m =e i7r ^i la ? a 'a m , m = 2,...,7V, (4.7) 
Eq. (|4.6p becomes a bilinear form in Fermi operators 



^7 M T A V- 



Jxyc^On 



4 



- ^(V -C 



) ( C L+1 



Cm+1 



(4.8) 



We assume periodic boundary conditions in Eq. 
omitting an unimportant boundary term in the limit 
TV — > oo for further calculations. To bring the Hamil- 
tonian to a diagonal form we perform first the Fourier 
transformation, 



1 H 

Trff E 



c — > p IKm r 

/tv 7 



(4.9) 



k = 27T//7V, / = -7V/2, . . . ,7V/2 - 1 (we assume TV is 
even), and then the Bogolyubov transformation, 



(3 K = i sin g K c K + cos g K J_ h 
J\ sin k 

tan(2g K ) = 



2 Jxy — Ji cos k 
The Hamiltonian (14.81) becomes 



(4.10) 



H = -^J H A + ^A K ^t/3 K -0, 



Jxy + ^Ji) - MxyJi cos2 ~. (4.11) 



Since Eq. (|4.11[) is the Hamiltonian which acts in the sub- 
space to which the ground state of the spin model (|4.ip . 
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(|4.5[) . (|4.3p belongs we immediately find the required 
ground-state energy per site in the thermodynamic limit 



eo = -g^nA _a_ E ( 2! ) ) 

2 2JxyJ\ 



(JXY + \Jl) 



(4.12) 



Here E(z) is the complete elliptic integral of the second 
kind (see Sec. with the modulus z given by (|4.12[) . 

It is worth to recall some ground-state properties of 
the spin model defined through the Hamiltonians (|4.ip . 
(|4.3[) . The energy spectrum A K is gapped unless 



Ji — 2 Jxy — 2 



D 2 . 



(4.13) 



If the condition (|4.13p holds the spin model has a gapless 
excitation spectrum^. From Eq. f|4. 13[) we see that only 



the relation between the Ising interaction J\ and the effec- 
tive XY component of the Heisenberg interaction Jxy is 
relevant for adjusting to a critical point, whereas the lat- 
ter effective interaction Jxy can be altered by varying 
the Dzyaloshinskii-Moriya interaction D. Furthermore, 
the condition (|4.13|) yields z 2 = 1 in Eq. (|4.12j) thus im- 
plying peculiarities in the ground-state properties of the 
investigated bond alternating chain. 

Now we turn to a deformable spin-1/2 Heisenberg- 
Ising chain with Dzyaloshinskii-Moriya interaction and 
assume that the interspin interactions depend on the 
change of intersite distance as follows: Jh = J7h(1 — k^S) , 
D = V(l — k D S), Ji = Ji(l - kiS). We introduce the 
parameters /ch, &d and k\ to distinguish between the ef- 
fects of the distance change on different intersite inter- 
actions. Using Eq. (|4.12p we may obtain the magnetic 
contribution eo(5) to the variational enthalpy per site 
g(T = 0,p;5) 



eo(S) = -IjhH 1 - fc H^) 



2 V^h(! - fe H<5) 2 + £> 2 (1 - k B S) 2 + Jt(l - hS) 



47T 



E(z), 



8V# - k H Sy + V 2 (l - k D 6) 2 - hS) 



V^ 1 - kaS) 2 +V 2 {l - k D S) 2 + Jj(l - k x S) 



(4.14) 



It is worthy to notice that the expression (|4.14[) for the 
ground-state energy of the Heisenberg-Ising bond alter- 
nating chain is from the mathematical viewpoint similar 
to the expression p. 21) determining the ground-state en- 
ergy of the transverse Ising chain. To reduce the number 
of free parameters we will further assume fcjj = &d k\ 
[note that the choice ku = ku = k\ implies that z in 
Eq. (|4.14[) becomes (^-independent and thus, this choice 
leads just to a trivial case with x = a]. Moreover, we 
set Jn = 1, A = 1 and kn = kr> = 1, hence leaving 
only three free parameters which control the interaction 
strengths and the effect of intersite distance change on 
the interspin interactions, Ji, T> and k\ = k. Under these 
assumptions, the expression (|4.14[) simplifies to 



... 1-6 2VTTW(l-S)+J 1 (l-kS) 
e (5) = i -i i iE(z) 



-in 



8Vl+V 2 J l (l-S)(l-k5) 
[2^1+^(1 - S) + J I (l-kS)\ 



(4.15) 



The value of S(p) and the inverse compressibility x are 
given by Eqs. (|3.3|) and (|3.4|) . respectively, however, with 
eo(S) now given by Eq. (|4.15j) . 
In what follows we will discuss the effect of an ap- 



plied pressure p on elastic properties of the quantum 
Heisenberg-Ising bond alternating chain defined through 
the Hamiltonians (|3~T]) . (|3~2"]) . (|Q|) . It should be men- 
tioned here that the pressure may generally induce dif- 
ferent changes in the relevant interaction constants and 
hence, the pressure variations may be used for adjust- 
ing the interaction constants so as to achieve a spe- 
cific condition inherent to a quantum critical point. In 
Figs. [51 [HI and [7] we demonstrate the dependences of 
the inverse compressibility x on the pressure p along 
with the relevant changes of 8{p) at different values of 
Ji = 1.8, 2, 2.2, k = 0, 2, and V = 0, 0.5, which 
were obtained by solving Eqs. (|3.3p . Q3.4p with eo(<5) 
(|4.15p by subsequent numerical integration and differ- 
entiation. It is quite clear that peculiarities in various 
quantities may be expected if the minimum of the vari- 
ational enthalpy g(T = 0,p;8), which occurs at 5(p) 
given by Eqs. (|3.3p and (|4.15p . satisfies the relation 
Ji(l-k5(p)) = 2(1- <J(p))V '1 +V 2 , i.e., when 6(p) = S c , 



2V1 + V 2 - Jx 



(4.16) 



2VT+V 2 - kjj 
By employing Eq. (|4.16p one may subsequently calculate 
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FIG. 5: (Color online). Inverse compressibility x(p) for the 
spin- 1/2 Heisenberg-Ising chain (|4.1j) . (|4.2|) . (|4.3j) with Jr - 
1, kn = 1, A = 1, T> = 0, fa = 0, a = 2 and three different 
values of Ji = 1.8 (red), Ji = 2 (blue), Ji = 2.2 (green). 
Inset: Corresponding S(p) versus p dependences. 




FIG. 7: (Color online). Inverse compressibility yc(p) for the 
spin-1/2 Heisenberg-Ising chain (T4~Tj) . (|4~2|) . (|4~3|) with Jk = 
1, fc H = 1, A = 1, V = 0.5, fc D = 1, ki = 2, a = 2 and three 
different values of Ji = 1.8 (red), Ji = 2 (blue), Ji = 2.2 
(green). Inset: Corresponding 5(p) versus p dependences. 




FIG. 6: (Color online). Inverse compressibility >c(p) for the 
spin-1/2 Heisenberg-Ising chain P~Tj) . P~2")) , (TO)) with J H = 
1, kn = 1, A = 1, T> = 0, fa = 2, a = 2 and three different 
values of Ji = 1.8 (red), Ji = 2 (blue), Ji = 2.2 (green). 
Inset: Corresponding 5(p) versus p dependences. 



from Eq. 



the critical pressure 



l _ 2yi + v 2 + fcjj 

8 4tt 



a — . (4.17) 

2VT+V 2 - kJi 



As one can see from the insets in Figs. [5l HI and [71 the 
most robust changes in the monotonous 5(p) dependences 
can be observed in the neighborhood of p = p c . Note 
furthermore that the greater the difference between the 
constants fcn = &d and k\ is, the more robust are the 
respective deviations from a linearity in the S(p) versus 
p dependence. Because of likeness of Eq. (14.151) and Eq. 
the behavior in the vicinity of 5 C is similar to the 



one treated in Sec. IIII| see the discussion around Eqs. 
(TXa - (rXTD . 

To gain an insight into the overall behavior, Fig.[S]illus- 
trates the situation when the pressure- induced changes of 
the Hcisenberg interaction are greater than that of the 
Ising interaction on behalf of fcu > fci, while the reverse 
case is displayed in Fig. [5] when considering another pos- 
sible particular case with /ch < k\. Altogether, it could 
be concluded that the quantum spin system placed on the 
linear-elastic lattice is responsible for the deviations from 
Hooke's law, i.e., the extension (or contraction) of the lat- 
tice is not linearly proportional to the applied pressure 
not because of anharmonicity of the lattice but owing to 
a critical behavior of interacting quantum spins placed 
on its sites. 

Finally, let us turn our attention to the inverse com- 
pressibility x{p) shown in the main panels in Figs. 
El and [7J It can be directly observed from Figs. [5] and 
[5] that the inverse compressibility exhibits at sufficiently 
low pressures a more intriguing pressure dependence with 
a substantial decline of the inverse compressibility inher- 
ent to the disordered phase (the regime of weak Ising in- 
teraction), while it always monotonically increases with 
rising pressure towards its maximum value x = a in 
the long-range-ordered phase (the regime of strong Ising 
interaction). Moreover, the quantum critical point is 
shifted towards higher pressures with increasing a rel- 
ative strength of the Ising interaction in the particular 
case with kji = ku > ki, whereas the opposite trend (i.e., 
a shift towards lower pressures) can be observed in the 
particular case with fcn = kjj < ki. Obviously, the ef- 
fect of the introduced Dzyaloshinskii-Moriya interaction 
does not fundamentally affects the physical picture we 
drew except that the relevant shift of quantum critical 
point due to the same change of the Ising interaction be- 
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comes much more pronounced as it can be seen from a 
comparison of Figs. [5] and [7] The most interesting find- 
ing stemming from the present study is that the inverse 
compressibility k(j>) continuously diminishes while ap- 
proaching the critical pressure p c quite similarly as does 
the inverse compressibility of the transverse Ising chain, 
compare Figs. [SJ [H and[7]with Fig.@] This result is taken 
to mean that the elastic properties of the lattice give a 
clear evidence for a quantum critical behavior of the spin- 
1/2 Heisenberg- Ising bond alternating chain (|4.1[) , (|4.2p . 

631). 

V. CONCLUSIONS 

Dcformablc spin lattices are often discussed in theoret- 
ical studies of, e.g., magnetic^ or hydrogen-bonded ferro- 
/antiferroelectrio2£ compounds. However, in those stud- 
ies various approximations are used which may mask de- 
tailed effect of spin-lattice coupling. In the present work, 
we have followed rigorously how lattice elastic constant 
manifests peculiarities of the spin system in the ground 
state driven by varying of an external field or pressure. 
More precisely, we have discussed some ground-state 
properties of several dcformablc quantum spin chains fo- 
cusing on the spin- 1/2 Hcisenberg-Ising bond alternat- 
ing chain refined by the antisymmetric Dzyaloshinskii- 
Moriya interaction. In particular, we have followed how 
peculiarities in the quantum spin systems manifest them- 
selves in the elastic properties under simple assumption 



about the spin-lattice interaction. We have found that 
the observed elastic constant k reflects in a rather obvi- 
ous manner (shows a decline) the behavior of magnetic 
subsystem in the vicinity of a quantum critical point. Al- 
though our analysis refers to the T = case, for which 
the studied effects are the most pronounced, it may be 
generalized in a rather straightforward way for nonzero 
temperatures T > too. 

The magnetic-field dependence of the compressibility 
at low enough temperatures is amenable to experimental 
measurements as for instance shown in Ref. [26| reporting 
such measurements for the spin-chain material CuCCP, 
or in Ref. [2?] reporting such measurements for the two- 
dimensional frustrated antiferromagnet CS2CUCI4. From 
this perspective, it is worthwhile to remark that the elab- 
orated analysis may be extended also for the spin-1/2 
Heisenberg chain after introducing a Hartree-Fock ap- 
proximation for Jordan- Wigncr fcrmions. 
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